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Abstract. The neutrino oscillations in Earth matter introduce modulations in the 
supernova neutrino spectra. These modulations can be exploited to identify the 
presence of Earth effects on the spectra, which would enable us to put a limit on 
the value of the neutrino mixing angle ^13 and to identify whether the mass hierarchy 
is normal or inverted. We demonstrate how the Earth effects can be identified at a 
single detector without prior assumptions about the flavor-dependent source spectra, 
using the Fourier transform of the "inverse-energy" spectrum of the signal. We explore 
the factors affecting the efficiency of this method, and find that the energy resolution of 
the detector is the most crucial one. In particular, whereas water Cherenkov detectors 
may need a few ten thousand events to identify the Earth effects, a few thousand 
may be enough at scintillation detectors, which generically have a much better energy 
resolution. A successful identification of the Earth effects through this method can 
also provide Amg to a good accuracy. The relative strength of the detected Earth 
effects as a function of time provides a test for supernova models. 
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1. Introduction 

Our knowledge of neutrino masses and mixing parameters has been rapidly improving 
in the past few years. We know that the neutrino flavors u^, and i/^ mix among 
themselves, and we already have measured the mass squared differences and two of 
the three mixing angles to a good accuracy [l]-[7]. The most important unknowns are 
the value of the third mixing angle, whether the neutrino mass hierarchy is normal or 
inverted, and the magnitude of the leptonic CP violation. The neutrino spectra from a 
core collapse supernova (SN) can shed light on the first two of these unknowns. 

Neutrinos produced inside a SN core undergo oscillations on their way out through 
the mantle and envelope of the star, through the interstellar space, and possibly even 
through some part of the Earth before arriving at the detector. The spectra of these 
neutrinos carry information about the two mass squared differences and the z/g flavor 
component in the three mass eigenstates. Of course this information comes convoluted 
with the primary fluxes of the neutrinos produced inside the star, and the extraction 
of the oscillation parameters depends crucially on our understanding of these primary 
fluxes. 

The uncertainties in the calculations of the primary flux spectra remain large, so 
only some of the robust features of these spectra can be used with confidence to extract 
the mixing parameters. Even with this limitation, it has been argued [8, 9] that the 
observations of the i^e and Ug. spectra at the detectors on Earth may reveal the type of 
the neutrino mass hierarchy and imply a range for the third mixing angle. Significant 
modifications of neutrino spectra can take place if the neutrinos travel through the Earth 
matter before reaching the detector, and these Earth effects can provide perhaps the 
most concrete signatures for some of the neutrino mixing scenarios [10, 11, 12]. 

Any signature of neutrino oscillations depends on the flavor-dependent differences 
between the primary spectra. Recent studies [13, 14, 15, 16] reveal that these differences 
are much smaller than had been assumed before. This makes the identiflcation of Earth 
effects harder than that expected in the previous analyses, and one needs to reevaluate 
the potential of the detectors from the perspective of these "pessimistic" assumptions 
about the primary fluxes and spectra. 

The comparison of the neutrino signals observed at two detectors would be the most 
efficient way of detecting the Earth effects. However the detectors need to be sufficiently 
far apart so that the Earth effects are different for both of them, and sufficiently large 
to observe a statistically significant signal. Currently Supcr-Kamiokande is the only 
large detector that can measure the neutrino energies and can detect more than a few 
thousands of events from a galactic SN. 

Though two detectors, both capable of measuring the neutrino energies, is definitely 
the most desirable option, it is still possible to detect Earth effects without having to 
measure the energies of individual neutrinos. Indeed, it has recently been shown that for 
a galactic SN, the comparison of the number of Cherenkov photons detected at Super- 
Kamiokande and IceCube as a function of time may be able to identify the Earth effects 
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[17]. The location of the SN can be anywhere within a specific 70% region of the sky 
for the effects to be observable. 

The motivations to look for a way to identify the Earth effects at a single 
detector are manifold, the cancellation of systematic uncertainties being the major 
one. Moreover, in conjunction with the two-detector signature mentioned above, the 
fraction of the sky for a favorable location of the SN increases to about 85% when Super- 
Kamiokande is also employed as the "single" detector. For nearly 35% of the sky fraction 
for the SN location, these two ways of detecting the Earth effects act as confirmatory 
tests of each other. The ranges of primary neutrino fiux parameters that give an Earth 
effect signature through these two procedures also differs slightly, since the two-detector 
method involving IceCube relies on the measurement of integrated luminosity, whereas 
the one-detector method relies on the analysis of the spectral shape. In this sense, the 
two methods are complementary to each other. 

The naive method of fitting the observed neutrino signal for the primary fluxes 
and neutrino mixing parameters is inefficient for several reasons. The primary fluxes 
- their average energies, spectral widths as well as the magnitude of the total fluxes - 
are time dependent, and only a few rough features can be said to be known with any 
confldence. One may try to get rid of the time dependence by dividing the signal into 
several time bins, but this reduces the available statistics. Even within each time bin, 
one typically has to fit for as many as 5 parameters of the primary spectra, not counting 
the uncertainties in the neutrino mixing parameters. 

However, robust identification of Earth effects can be achieved by observing that the 
parameters that govern the oscillation frequency of the neutrinos inside the Earth are 
relatively well measured, and more importantly, completely independent of the primary 
neutrino spectra. Therefore, the Earth effects can be identified merely by identifying 
the presence of this oscillation frequency in the observed spectrum. The oscillation 
frequency does not change with time, thus obviating the need for several time bins. This 
oscillating component is expected to be a small addition to the otherwise approximately 
blackbody spectrum that forms the major component of the signal. In order to extract 
this small oscillation component, we propose a Fourier analysis of the spectrum, which 
can separate the signals of different frequency from their superposition. We perform 
a numerical simulation to explore the efficiency of this method at a water Cherenkov 
detector like Super-Kamiokande, and a large scintillation detector like the one proposed 
in Ref. [18], which has a much better energy resolution. 

Since the Earth effects can be observed only for certain neutrino mixing scenarios, 
the mere identification of these effects can already tell us whether, for example, the 
neutrino mass hierarchy is normal. This is a daunting task even at the future long- 
baseline experiments. In addition, the magnitude of the Earth effects as a function 
of time gives us important information about the primary neutrino spectra that is 
extremely useful for understanding the SN dynamics. 

This paper is organized as follows. In Sec. 2, we describe the modifications of 
the primary neutrino spectra by the Earth matter effects. In Sec. 3, we introduce the 
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Fourier transform for extracting the oscillation frequency of neutrinos inside the Earth, 
and propose a procedure for identifying the Earth effects without having to make any 
assumptions about the primary neutrino spectra. In Sec. 4, we explore this method at 
scintillation and water Cherenkov detectors through a Monte Carlo simulation. Sec. 5 
concludes. 



2. Neutrino mixing parameters and Earth effects on the Ue spectrum 

Neutrino oscillations are now firmly established by measurements of solar and atmo- 
spheric neutrinos and the KamLAND [19] and K2K [20] long-baseline experiments. The 
weak interaction eigenstates Vg, and v-j- are non-trivial superpositions of three mass 
eigenstates i^i, V2 and 1/3, 

^a^U^iUi, a e {e,//,T}, i e {1,2,3} , (1) 

where U is the leptonic mixing matrix. It can be written in the canonical form 

U — -R23(^23)-Rl3(^13)-Rl2(^12) , (2) 

where RijiOij) corresponds to the rotation in the i-j plane through an angle 6ij. We 
have neglected the CP violating effects, which are irrelevant for SN neutrinos. Prom a 
global 3-flavor analysis of all the available data, one finds the Scr ranges for the mass 
squared differences Am?- = mf — rrij and mixing angles as summarized in Table 1. 



Table 1. Neutrino mixing parameters from a global analysis of all experiments 
(3(7 ranges) [7]. 



Observation 


Mixing angle 


Am^ [meV^] 




Sun, KamLAND 


012 = 32°- 


-42° 


Am^x — rn\ — m\ ~ 


54-190 


Atmosphere, K2K 


023 = 34°- 


-60° 


lAmi^l = \ml~ml 


= 1500 3900 


CHOOZ 


^13 < 13° 




Arn^i = TO§ — nil ~ 


Aml2 



A SN core is essentially a neutrino blackbody source, but small flavor-dependent 
differences of the fluxes and spectra remain. Since these differences are very small 
between and P^, we represent both these species by Px- We denote the fluxes of Pe 
and Pj; at Earth that would be observable in the absence of oscillations by F° and F°, 
respectively. Por both Fg and F^, we assume a distribution of the form [15] 

E 



^^^^~Eo r{l + a) \Eo) 



(3) 



where Eq is the average energy, a a parameter that typically takes on values 2.5-5 
depending on the flavor and the phase of neutrino emission, and $o the overall flux at 
the detector. The values of the total flux $0 and the spectral parameters a and Eq are 
different for z/g and u^- We represent these by the appropriate subscripts. Two of the 
most important robust features of the primary spectra are 

(1) the energy hierarchy of the neutrino species: {E^^) < {Ep^) < (Ep^), 
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(2) "pinching" of the spectra: ctj > 2 for all species. 
Apart from a heavy water detector like SNO, the neutrino detectors are sensitive mainly 
to z/g in the SN energy range. We shall therefore concentrate only on the z/g spectrum 
in this paper. In the presence of oscillations a z/g detector actually observes the flux 

Fi>{E) = p''{E)El{E) + [l - p^iE)] F^iE) , (4) 

where p^{E) is the Pg survival probabihty after propagation through the SN mantle 
and perhaps part of the Earth before reaching the detector. The bulk of the Pg ^-re 
observed through the inverse beta decay reaction Pgp — > ne"*". The cross section a of 
this reaction is proportional to Ep , making the spectrum of neutrinos observed at the 
detector N{EpJ oc aF^ oc ElF^. 

A significant modification of the survival probability due to the propagation through 
the Earth appears only for those combinations of neutrino mixing parameters shown in 
Table 2. The Earth matter effects depend strongly on two parameters, the sign of Am|2 
and the value of sin^ ^13 [8, 10]. The "normal hierarchy" corresponds to m-i <m2 < rus, 
i.e. Am32 > 0, whereas the "inverted hierarchy" corresponds to < mi < m2, i.e. 

< 0- Note that the presence or absence of the Earth effects discriminates between 
values of sin^ ^13 less or greater than 10^^, i.e. |^i3| less or larger than about 1.8°. Thus, 
the Earth effects are sensitive to values of 16*131 that are much smaller than the current 
limit. 



Table 2. The Earth effects appear for the indicated flavors in a SN signal. 



13-Mixing 


Normal Hierarchy 


Inverted Hierarchy 


sin^ < 10-3 


and 


Vg and Vg 


sin^ 6I13 > 10-3 







Let us consider those scenarios where the mass hierarchy and the value of ^^13 are 
such that the Earth effects appear for Pg. In such cases the Pg survival probability p^{E) 
is given by 

^ cos^ 012 - sin 2^® sin(2^® - 2^12) sin^ (^12.5 j , (5) 

where the energy dependence of all quantities will always be implicit. Here 9% is 
the mixing angle between Pg and P2 in Earth matter while Am^ is the mass squared 
difference between the two anti-neutrino mass eigenstates Pi and P2 in units of lO^^eV^, 
L is the distance traveled through the Earth in units of 1000 km, and E is the neutrino 
energy in MeV. We have assumed a constant matter density inside the Earth, which is a 
good approximation for L < 10000 km, i.e. as long as the neutrinos do not pass through 
the core of the Earth. 

The energy dependence of p^ introduces modulations in the energy spectrum of Pg. 
These modulations may be observed in the form of local peaks and valleys in Fig. 1(a), 
which shows the spectrum of the event rate, aF^, at a detector as a function of the 
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neutrino energy. Fig. 1(b) shows the same neutrino signal as a function of the "inverse- 
energy" parameter, defined as 

y=12.5/E . (6) 

Whereas the distance between the peaks of the modulation increases with energy in 
the energy spectrum, the peaks in the inverse-energy spectrum are nearly equispaced 
and hence have a single dominating frequency. This makes it easier to distinguish these 
modulations from random background fluctuations that have no fixed pattern. 




Figure 1. The energy spectrum (a) and the inverse-energy spectrum (b) of aPf . 
The fluxes are normalized such that the area under each curve is unity. For all the 

examples in this paper, we use the primary neutrino flux parameters ap^ = ap^ = 3.0, 
{EpJ = 15 McV, {£",7,,) = 18 McV, = 0.8, which arc realistic for the fluxes 

during the cooling phase. For the mixing parameters, we use AttIq = 6 (in 10^^ eV^) 
and sin^(2^0) = 0.9. The distance travelled through the Earth is L = 6 (in 1000 km) 
unless otherwise specified. 

The equidistant peaks in the modulation of the inverse-energy spectrum are a 
necessary feature of the Earth effects. Indeed, the net Aux at the detector may 
be written using (4) and (5) in the form 

= sin^ OuF^ + cos^ ^12^° + AF^i® sin' (Am|Ly) , (7) 

where AF^ = {Fg — F^) depends only on the primary neutrino spectra, whereas 
= —sin 2(9^ sin(2^^ — 26*12) depends only on the mixing parameters and is 
independent of the primary spectra. The last term in (7) is the Earth oscillation term 
that contains a frequency = 2Am0L in y, with the coefficient AF^A^ being a 
comparatively slowly varying function of y. The first two terms in (7) are also slowly 
varying functions of y, and hence contain frequencies in y that are much smaller than 
k^. The dominating frequency A;^ is the one that appears in the modulation of the 
inverse-energy spectrum in Fig. 1(b). 

The frequency k^ is completely independent of the primary neutrino spectra, and 
indeed can be determined to a good accuracy from the knowledge of the solar oscillation 
parameters, the Earth matter density, and the direction of the SN. If this frequency 
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component is isolated from the inverse-energy spectrum of Pg, the Earth effects would 
be identified. In the next section, we shall show how this can be achieved through the 
Fourier transform of the inverse-energy spectrum. 

3. Identifying the Earth matter effects 

3.1. Fourier transform of the inverse-energy spectrum 

Taking the Fourier transform of a function is the standard way of extracting components 
of different frequencies present in the function. The Fourier transform of a function f{y) 
is 

/CXD 
fiy)e'''dy (8) 
-oo 

while the "power spectrum" Gf{k) = gives the strength of the frequency k 

present in f{y). 

In Fig. 2(a), we show the power spectrum G„F{k) of the y-spectrum aFP . The peak 
at /c ~ 2AmQL = 72 corresponds to the oscillations in Earth matter. The large peak 
at low values of k, which has the value of unity at /c = 0, is the dominant contribution 
due to the first two terms in (7). As may be observed, the Earth effect peak is cleanly 
isolated from the dominant contribution. The figure also allows us to put a lower bound 
on the value of k^ for which the Earth effects will be detectable: if Am^L < 20, the 
Earth effect peak will be lost in the dominant low frequency peak. Physically speaking, 
this implies that the neutrinos have to travel a minimum distance through the Earth for 
the Earth effects to be detectable. So this is not a limitation of this particular method, 
but a general limitation on the identification of Earth effects. 




20 40 60 80 100 20 40 60 80 100 



Figure 2. The power spectrum GcF{k) of aFf (a) and the same power spectrum 
with a sharp energy threshold of 10 MeV (b). 

The Earth effect peak has a finite width due to the finite ^/-dependence of the 
coefficient of the oscillating term AF^^?®, and of Am^. The magnitude of the width 
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is a measure of the ^/-dependence of these quantities, whereas the area under the peak 
gives the total contribution of the oscillating Earth effect term. 

In Fig. 3(a), we show the envelope of the Earth effect term, = aAF^Af^, for 
the typical neutrino spectra parameters during the cooling phase, normalized such that 
/ aF^dy = 1. The power spectrum of the Fourier transform Gxijy) of this envelope is 
shown in Fig. 3(b). Since is the coefficient of the Earth matter oscillations term 
siv?{k(^y) and /c® is nearly constant over all the energy range, the width of the Earth 
effect peak in Fig. 2(a) is due mainly to the width of Gt(^)- We parameterize this width 
by defined as the value of k where GT{k) is half its maximum value. This would 
correspond to the half width at half maximum of the /c© fa 72 peak in Fig. 2(a). 




.2 . 4 . 6 . 8 1 1 .2 1 . 4 10 20 30 40 50 



Figure 3. The Earth term = aAF^A^ as a function of y, and its power spectrum 
Grik). 

In Fig. 3(b), we observe w ~ 7. The value of w depends on the primary fiuxes, but 
it is in the range of 3-10 for almost all the allowed parameter range. The exact value of 
w is not crucial for this analysis, however this estimation is useful in two ways. It allows 
us to get rid of spurious peaks by applying a selection criterion of a minimum width of 
the peak, and it also gives us the maximum range of /c- values around the mean for 
which the Earth effect term contributes. 

The neutrino detectors have an energy threshold of about 5-10 MeV, which 
corresponds toy = 1.25-2.5. Although the number of events below this energy threshold 
is expected to be very small (Sec Fig. 1), the presence of a sharp threshold may introduce 
high frequency components in the y-spectrum. Although the threshold is not sharp in 
real life, we show an extreme example of the effect of a sharp threshold in Fig. 2(b), 
taking a threshold of 10 MeV. Since the Fourier transform of a step function is a series of 
equispaced frequencies, we observe that these high equispaced frequencies are superposed 
on top of the actual Gapik). The spacing between these frequencies is much smaller than 
the width of the Earth effect peak, so the observation of the peak and the area under 
the peak are not much affected, even in this extreme situation. In reality, the energy 
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threshold is expected to be much smoother and should not affect the identification of the 
peak and the measurement of its strength. In all the numerical simulations henceforth, 
we introduce a sharp energy threshold of 5 MeV in order to take care of the threshold 
effects. 

3.2. The background due to statistical fluctuations 

Since the spectrum is observed as a discrete set of neutrinos with individual energies 
(and hence, individual y values), we have to approximate g{k) by the discrete sum 

9{k) = ^ E e^'^ > (9) 

V » events 

where is the total number of events. The statistical fluctuations in the spectrum 
can contribute to frequency components that are comparable to A;®. This acts as a 
background for the actual Earth effect peak. 

The magnitude of this background can be estimated by assuming that the 
distribution of ky (modulo 27r) is completely random. This is definitely a good 
approximation for large values of k independent of the details of the distribution of y. 
In this approximation, the sum in (9) above represents a two dimensional random walk 
with unit step. This implies that the value of GN{k) obeys an exponential distribution 
P{G) = e^'^, which has a mean fi = 1.0 and variance = 1.0 independent of the 
number of events N. The background due to the statistical fluctuations is thus nearly 
1.0 on average, and the "signal" of the peak should rise above this average level in order 
to be identified. 

This is illustrated in Fig 4(a), where the distance travelled by the neutrinos through 
the Earth is L = 0. Note that the low frequency background is dominant for k < 40. In 
Fig. 4(b), we show the same neutrino signal, but now with L — 6. The Earth effect peak 
stands out above the background. Though the level of the background is not affected by 
the number of events, the magnitude of the signal grows proportional to N. Clearly, the 
Earth effect peak would grow above the background level when the number of events is 
sufficiently high. 

Although the peak identification may be possible "by eye" in many cases, as 
in Fig 4(a), the standard particle physics method of resonance identification may be 
employed for the task. Firstly, we apply a lower /c-cut to cut away all k < kcut, where 
kcut is the value of k where the dominant peak falls to the background level. In Fig 4(a), 
for example, kcut ~ 40. Although we have made a rough estimate of the background 
above, the actual background may be calculated by averaging the spectrum Gat (A;) away 
from the value fc®, i.e. in the region where we do not expect the Earth effect peak. 

We know from Fig. 2 that the half width at half maximum of the peak in k should 
be w ~ 3-10. We then expect that most of the Earth effects contribute to values of k 
within a range of W around the peak, where W may be taken to be 4w, for example. 
We can then choose any interval in k with a range of W, and integrate the power 
spectrum within this range. Applying the Central Limit Theorem to the background 
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(a) 

L=0 km 





(b) 

L=6000 km 
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Figure 4. Background due to statistical fluctuations when the Earth efi^ects are absent 
(a) and the signal peak with the Earth effects (b). 



distribution P{G), we can estimate the integral of the background within this interval 
to be 5 = i^bW ± aeVW. For i^b = 1-0, aB = 1.0 and W = 20, we get 5 20 ± 4.5. 
This implies that /ly GAr(/c) dk > 35, for example, would correspond to a more than 3a 
detection of a positive signal. Depending on the actual data, the value of W may be 
optimized. The procedure of measuring the area under the peak is more efficient than 
just measuring the peak height since it can weed out the spurious high peaks that do 
not have the minimum width dictated by Gt(^)- 

The number of events needed in order to see the signal above the background 
depends strongly on the primary neutrino spectra, since the coefficient of the Earth 
effect term depends on AF". In Fig. 5(a), we show the Earth effect term sin^(/c®|/) 
in comparison with the term aF^, normahzed such that / crF| = 1. In Fig. 5(b) we 
show the normalized term aF^ for comparison. The area under this curve is unity. 
The area A under the curve in Fig. 5(a) may then be taken to be a measure of the 
net contribution per event of the Earth effect term to the y-spectrum. For the primary 
spectrum parameters used here, A ~ 0.03 per event. 

We show in Fig. 6 the dependence of A on the parameters of the primary spectra. 
We keep the parameters ap^ = 3 and (Ep^) = 15 McV fixed, and show ^ as a function of 
<I>P^/$P^ and (Ep^) for two typical values of a^^^. During the accretion phase, ap^ > ap^, 
which is the situation depicted in Fig. 6(a). During the cooling phase, the a values 
for both spectra are very similar, which is depicted in Fig. 6(b). The stars roughly 
correspond to the parameter values obtained in the SN simulations [16]. It may be 
observed that for these parameter values, A <^ 0.01 during the accretion phase and 
A ~ 0.03 during the cooling phase, so that the Earth effects during the accretion phase 
are expected to be much smaller. The relative magnitude of these effects thus also 
provides a test for the SN simulations. 
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Figure 5. The oscillating Earth effect term, with the cross section of the inverse beta 
decay factored in (a). It is normalized such that / aF^ = 1. The normalized term aF^ 
is shown in (b) for comparison. Note that the scale in (b) is 10 times the scale in (a). 



(a) (b) 




15 16 17 18 19 20 15 16 17 18 19 20 




Figure 6. The dependence of A, the area under the curve in Fig. 5(a), on the primary 
neutrino flux parameters. The contour lines have values starting from 0.05 at the 

lightest end, with decrements of 0.01 going towards the darker regions. The panel (a) 
uses the parameters {Ep^) = 15 MeV, = 4, = 3, panel (b) uses {Ep^) = 15 
MeV, ap^ = ap^ = 3. 

We can estimate the number of events required for a peak identification as follows. 
In order to get a 3(7 identification, we need the "signal" contribution S such that 
S > 3^(7| + where erf = S. With cr^ = 4.5, this corresponds approximately to 
S > 20. In order to reach this strength S with an average contribution oi A~ 0.03 per 
event for example, we need at least S/A^ 700 events. This is a very rough, and indeed 
an optimistic, estimate of the minimum number of events needed in order to detect 
Earth effects unambiguously. In reality, the energy resolution of the detectors tend to 
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smear out the oscillations and decrease the magnitude of their strength, thus increasing 
the required number of events. We shall study the effects of the energy resolution of 
the detectors in Sec. 4 and get a realistic estimate of the number of events needed to 
identify the Earth effects. The relative strength of the Earth effects for different primary 
spectra can still be read off from Fig. 6. 

Note that the above procedure can be employed without any prior knowledge of k^. 
Actually, if the value of Am^ and L is known, we already know the value of to look 
for. This helps in getting rid of any background due to spurious peaks. On the other 
hand, if this peak is identified unambiguously, the value of k^ can help in improving the 
accuracy of the measurement of AtHq. We shall study this in the next subsection. 

3.3. Determination of Am^ 

The current 3a range of the solar mass squared difference Am| = (5.4-19) x 10~^ eV^ 
is obtained mainly through the combination of the limits from Supcr-Kamiokande and 
KamLAND. Although this range is expected to narrow significantly with the future 
KamLAND data, it is worthwhile to note that the Fourier analysis of the SN neutrino 
spectra can also determine this value to an accuracy of a few percent. 

Once the Earth effect peak is identified, the value of A;® gives the value of Am^ 
since the value of L should be well known once the SN direction is estabhshed. The 
error in the measurement of the position of the peak may be roughly estimated hy w/k^. 
Since we expect w ^ 3-10, and k^ may be as high as 2Am|(max) L(max) 400, even 
this conservatively estimated error may be only a few percent. As long as A;© > 40, 
which is the minimum value of kq^ for the Earth effects to be detectable, the error due 
to determination of the peak position is less than 25%. 

The value of Am| is related to AtUq by 



where V is the magnitude of the matter potential inside the Earth. In Fig. 7, we show 
the y-dependence of Am|/Amg with various values of solar parameters. Since the y- 
spectrum is significant only for y > 0.2, the deviation of this ratio from unity at ~ 0.2 
may be taken as a conservative estimate of the error on Am^ from the SN spectral 
analysis. The figure shows that for y > 0.2, the values of Ami, and Am^ differ by less 
than 25%, and this difference decreases with increasing Am^ values. Therefore, it may 
be safely assumed that effectively, we have Am© ~ AtUq to within 25%. 

Throughout our analysis, we have assumed a constant matter density inside the 
mantle of the Earth. Actually the density may vary by as much as 30% along the 
neutrino trajectory in the extreme case where the trajectory is nearly tangent to the 
core. Since V (x p, this density variation contributes to the energy dependence of Am© 
through the term involving EV in (10), and smears the Earth effect peak. However, 
since the variation in by a factor of five is already taken care of in the above error 
estimation, the additional smearing due to the Earth density variation may be safely 
neglected. 



Ami = 



Am| sin^ 20q + (cos 20q + 2 VE/Amlf 



a 1/2 



(10) 
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Figure 7. The ratio Am^/Amg as a function of y for different solar parameters. 

6 X 10-5 eV2, sin2(20i2) = 0.7, 
2(2012) = 0.9, 

Dashed line: Am| = 18 x IQ-^ eV^, sin2(26ii2) = 0.7, 



Dot-dashed line: AtUq 

Dotted line: Am| = 6 x 10"^ eV^, sin2(26'i2) 



SoUd line: Atoq 



18 X 10-5 eV2, sin2(26ii2) = 0.9. 



Combining the error due to the peak determination and the error due to the 
difference in Am^ and Am^, even conservatively the value of Am^ would be known 
to within 35%. KamLAND may have already pinned this value down to a much better 
precision, however since no other planned experiment would be able to reach a 35% 
precision on AtUq in near future, this could be the first confirmation of the Am^ value 
measured at KamLAND. At high values of Am^ and L, the error due to the peak width 
as well as the error due to the Am| — Am| difference decreases and the accuracy in 
the AtUq measurement may become competitive with KamLAND. For example, with 
AtUq = 15 X 10~^ eV^ and L — 3000 km, the error in Am| from the SN spectral analysis 
would be less than 15%, comparable with the accuracy expected at KamLAND [21]. 

4. The effect of finite energy resolution 

The neutrino signal measured at the detector through the inverse beta decay process 
PeP — > ne'^, which is the dominant reaction in both the water Cherenkov as well as the 
scintillation detector, can be written in the form 

N,{E,) oc / dE', n{E,,E',)£iE',) J dE d^os 9^^^^^ F,{E) ,(11) 

where E is the energy of the incoming antineutrino, E'^ is the true energy and E^ is the 
measured energy of the outgoing positron. The differential cross section of the inverse 

beta decay process is [d'^a{E'^, E, cos 9)/ {dE'^ rfcos^)], where 9 is the angle between the 
antineutrino and the positron. We denote by S{E'^) the efficiency of positron detection, 
and TZ{Ee, E'^) is the energy resolution function. We neglect the elastic u-e scattering 
interactions, which account for less than a few percent of the total number of events, 
and can even be weeded out through their strongly forward peaked angular distribution. 
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We take S{E'^) to be uniform over the whole energy range except for a sharp lower 
threshold of 5 MeV. The resolution function TZ{Ee, E'^) is a Gaussian distribution with 
mean E'^ and the standard deviation ctd given by the energy resolution of the detector 
D, where D is SC for a scintillation detector and SK for a Cherenkov one. We take 
the inverse beta decay process to be completely elastic, so that the differential cross 
section involves a S function nonvanishing only &t E'^ ^ E /[I + {1 — cos 9)E/m^, where 
rUp is the proton mass. The angular distribution is approximated to be isotropic, and 
the cross section proportional to E'^. This corresponds to the differential cross section 
in [22] in the limit of degenerate neutron and proton masses, vanishing electron mass, 
and isotropic angular distribution. These approximations made here retain the essential 
features of the observed spectra. We finally normalize N^{Ee) to the number of events 
using /Are(^e) dEe = N. 

The outgoing positrons are detected in the water Cherenkov detectors through 
the Cherenkov photons that they radiate. In the scintillation detectors, the positrons 
are detected through photons produced in the scintillation material. Since a larger 
number of photons can be produced in a scintillation detector, these have typically a 
much better energy resolution than the water Cherenkov detectors. Indeed, whereas at 
Super-Kamiokande we have iTsk ~ 1-5 MeV (£^/10 MeV)^/^ [23], the energy resolution 
of a scintillation detector may be as good as ctsc ~ 0.22 MeV {E/IQ MeV)^/^ [18], 
which is more than a factor of 5 better. 

Since the identification of Earth effects as described in this paper relies on the 
detection of the modulations (see Fig. 1), and the finite energy resolution tends to smear 
them out, it is clear that the energy resolution plays a crucial role in the efficiency of 
detecting Earth effects. To illustrate this, we show in Fig. 8 the same energy spectrum of 
the signals as observed at a scintillation and a Cherenkov detector. One can clearly see 
how the poor energy resolution of the Cherenkov detector smears out the modulations. 
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Figure 8. The energy spectrum as observed at a scintillation detector (SC) and a 
Cherenkov detector (SK). All parameters except the energy resolution are identical. 
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The difference between Figs. 8(a) and (b) is reflected in the number of events 
required for the signal to rise above the background at the two detectors. A numerical 
simulation that generates inverse beta decay events in each of the detectors illustrates 
the comparative efficiency of the Fourier analyses at these two detectors. In Fig. 9, we 
show the Fourier transforms of the spectra at these two detectors with 2000 events each. 
Whereas the peak can be identified even by eye at the scintillation detector, the Fourier 
power spectrum at the Cherenkov detector is indistinguishable from background. 




25 50 75 100 125 150 25 50 75 100 125 150 

Figure 9. The Fourier transforms of the simulated y-spectra in scintillation and 
Cherenkov detectors with 2000 events each. 



In Fig 9(a), the area under the peak is around 80. Since the area under the peak 
grows nearly linearly with the number of events, this indicates that at a scintillation 

detector, the signal starts becoming visible at around N 1000. At the proposed 50 
kiloton scintillation detector LENA [18] for example, one expects about 13,000 events 
from a SN at 10 kpc. This is an order of magnitude more than the statistics required 
for identifying the Earth effects. Indeed, a much smaller scintillation detector may be 
sufficient for this purpose. 

At a water Cherenkov detector, the energy resolution is poor compared to a 
scintillation one. Indeed, the energy resolution is of nearly the same size as the 
wavelength of the Earth effect modulations for ~ 50. It is therefore much harder to 
identify the modulations. In Fig. 10(a), we show the Fourier transform of the y-spectra 
generated at a Cherenkov detector with 100,000 events. We expect that we need more 
than 60,000 events to be able to identify the peak unambiguously. This indicates that 
Supcr-Kamiokandc may be too small for detecting the Earth effects by itself in this 
parameter range. However the proposed Hyper-Kamiokande would have the required 
size. 

If the modulation wavelength is larger than the energy resolution, which is the 
case for low A;®, the peak is easier to detect. In addition, if the differences in the 
primary spectra of Pg and are larger than those taken in this paper till now, the peak 
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Figure 10. The Fourier transforms of the simulated ?y-spectra in a Cherenkov detector, 
(a) uses 100,000 events and the typical neutrino spectral parameters, while (b) uses 
more extreme "favorable" parameters, {Ep^) = 21 MeV and L = 4000 km. In the 
latter case, only 5,000 events are enough for the identification of Earth effects. 



identification can be achieved with a much smaller number of events. In Fig. 10(b), we 
show the power spectrum with the average energies of the primary spectra differing by 
40%, and the value of L decreased to 4000 km. In such a "favorable" range, even a 
few thousand events would be enough at a water Cherenkov detector, so that Supcr- 
Kamiokande has a chance of performing a successful Earth effect detection. However 
in the light of the latest calculations of SN spectra, such a large energy difference is 
basically ruled out. 

Note that even for a scintillation detector, the energy resolution becomes 
comparable to the modulation wavelength for large values, so the number of events 
required will increase with k^. At values of k^ > 200, the modulations get significantly 
washed out even at a scintillation detector. 

5. Conclusions 

The modulations introduced in the neutrino spectrum by the Earth matter effects 
provide a way of detecting the presence of these effects without prior assumptions about 
the flavor-dependent source spectra. We have demonstrated that a Fourier analysis of 

the inverse-energy spectrum of the z/g signal may reveal a peak corresponding to the 
neutrino oscillation frequency in Earth. The position of this peak is insensitive to the 
primary spectra, and depends only on Am^ and the distance travelled through the 
Earth. We study the feasibility of the identification of this peak. 

The number of events required for an unambiguous identification depends crucially 
on the energy resolution of the detector. The task is certainly feasible at a large water 
Cherenkov detector like Hyper-Kamiokande that can detect nearly 10^ events from 
a galactic SN. However, scintillation detectors generically have a much better energy 
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resolution than the Cherenkov ones, and typically a scintillation detector can detect the 
Earth effects with only about 2000 events. This indicates that a scintillation detector 
large enough to detect a few thousand events from a galactic SN is highly desirable. 

The identification of the Earth effects severely restricts the neutrino mixing 
parameter space, since the effects are present only with certain combinations of the 
neutrino mass hierarchy and the mixing angle ^13. In particular, if sin^^i3 is measured 
at a laboratory experiment to be greater than 10~^, then the Earth effects on the Pg 
spectrum imply the normal mass hierarchy. However if the Earth effects are not detected, 
it does not rule out any neutrino mixing parameters, owing to the current uncertainties 
in the primary fluxes. 

A galactic SN is a rare event, expected to occur only a few times per century. 
However this time scale should be compared with that of those laboratory experiments 
which would be sensitive to the mass hierarchy and values of sin^ ^^13 as low as 10^^. 
Determination of these two quantities is a difficult challenge even for the experiments 
that may be running in the next few decades [24, 25, 26]. If a SN is observed within 
the next decade, the information gained may even be useful in fine-tuning the design 
parameters of these experiments. 

The Fourier analysis described in this paper also yields the value of Atoq to a good 
accuracy if the matter effects are identified. Though KamLAND should pin down Am^ 
to within a few percent in the next few years, a galactic SN will provide a completely 
independent confirmatory test of this value. For large values of Ami, and L, the accuracy 
in Am|, through the SN spectral analysis becomes comparable with KamLAND. 

Models of SNe generically predict time variations in the neutrino flux spectra and 
hence in the magnitude of Earth effects. The relative strength of the detected Earth 
effects as a function of time thus provides a test for these models. 

In order to reap the benefits of the Earth effects, the detector needs to be shadowed 
by the Earth from the SN. The probability of such an occurrence is increased by having 
more than one detector spaced far apart from each other. From the point of view of SN 
neutrinos, several "small" (a few kiloton) scintillation detectors at different locations 
are far more useful than a single large one. 
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